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Conventional finite-temperature perturbation theory in 



which propagators have poles at k 



is shown to break 
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down at the two-loop level for self-interacting scalar fields. 
The breakdown is avoided by using free thermal propaga- 
tors that have poles at the same energy as the exact thermal 
propagator. This quasiparticle energy £(k) is temperature- 
dependent, complex, and gauge invariant. An operator theory 
containing two self-adjoint scalar fields is presented in which 
all temperature dependence is incorporated into the Hamil- 
tonian. No thermal traces are required to compute thermal 
Green functions. Choosing the spectrum of the unperturbed 
part of the Hamiltonian to contain the exact quasiparticle en- 
ergy £ (fe) produces a resummed perturbation theory that has 
the correct poles and branch cuts. The location of the poles 
and cuts is explained directly in terms of the spectrum of the 
Hamiltonian. 
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I. INTRODUCTION 

In high temperature QCD the self-energies of both 
quarks and gluons change the particle dispersion rela- 
tions from their zero-temperature form. The modified 
dispersion relations indicate that the naturally propagat- 
ing modes are collective excitations, i.e. quasiparticles. 
The locations of the poles in the propagators are gauge- 
fixing invariant [EJ . Braaten and Pisarski showed that for 
a consistent perturbation theory the propagators should 
not be free but should include part of the self-energies Q . 
With this hard-thermal-loop approximation to the quasi- 
particle effects it is possible to compute gauge invariant 
rates in which the corrections are controllably small, at 
least in certain kinematic regimes. A completely different 
test of the HTL quasiparticle dispersion relations comes 
from the excellent agreement with lattice data that is 
obtained by using an ideal gas of quark and gluon quasi- 
particles. When running coupling constants are incorpo- 
rated, the agreement holds even very close to the critical 
temperature of QCD ||. 

Quasiparticles are ubiquitous in finite-temperature cal- 
culations, and yet they normally play no role in the op- 
erator structure of the theory. This situation contrasts 
markedly with that of conventional particles at T=0. In 
particular, zero-temperature field theory has the follow- 
ing simple properties: 



• Poles in the propagator occur at the eigenvalues of 
the single particle state vectors. 

• Branch cuts in the propagator occur at the eigen- 
values of the multiparticle state vectors. 

This paper will describe a formalism in which these fea- 
tures also hold for finite-temperature quasiparticles. The 
present discussion will only treat scalar field theory. 

Sections II and III summarize the results of a recent 
paper on branch cuts at finite temperature which 
shows that the location of the cuts is determined by the 
quasiparticle energies. Section IV is based on a new for- 
mulation of finite-temperature field theory |^] in which all 
the Bose-Einstein temperature dependence is contained 
in the Lagrangian. Green functions are computed as vac- 
uum matrix elements rather than thermal traces. Sec- 
tions V and VI show how to split the Lagrangian into an 
unperturbed part that describes free quasiparticles and 
a remainder that determines the interactions. Several 
illustrations are worked out. The developments in Sec. 
IV and subsequently may be read independently of the 
earlier sections. 



II. BREAKDOWN OF PERTURBATION 
THEORY IN MASSIVE (f> 3 FIELD THEORY 

Braaten and Pisarski showed that finite-temperature 
perturbation theory breaks down in field theories with 
massless particles or those in which the mass is much 
smaller than the temperature. Resummation of hard 
thermal loops is then necessary for Green functions eval- 
uated at low energy and momentum. 

This section will discuss a more primitive breakdown 
of perturbation theory that occurs even without soft ex- 
ternal momenta and without high temperatures. The 
discussion will treat scalar field theory having cubic self- 
interactions, 
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• The Hamiltonian operator has eigenstates 
linear vector space. 
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with counterterms omitted. Perturbation theory is con- 
ventionally defined by choosing free thermal propagators 
that have poles at the zero-temperature mass m. This 
produces two-particle branch points in the self-energy at 
k 2 = Am 2 and three-particle cuts at k 2 — 9m 2 . Simple 
calculations show that at finite temperature perturbation 
theory fails in the vicinity of these points. The reason is 
that they are not the true branch points of the full theory. 
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A. Zero-temperature example 

A simple zero-temperature example will illustrate how 
higher order corrections can shift the location of branch 
cuts. Even though m is the physical mass for the theory, 
suppose that one performs perturbative calculations us- 
ing a free propagator A(fc) = l/[fc 2 — where mo is 
some different mass. It is important that mo is finite and 
not the bare mass. For definitcness choose m < m. The 
one-loop self-energy shown in Fig. 1 has a branch cut for 
k 2 > Am 2 ,. The discontinuity across the branch cut is 
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y^9(k 2 -4m 2 ). (2) 



Since m has no physical significance, this is clearly not 
a branch cut of the full theory. 



FIG. 1. One-loop self-energy. 



The indication that k 2 



Am 2 , is not a branch point 



of the full theory comes from the two-loop contribution. 
The full propagator is D'(k) = l/[k 2 — m 2 — II(fc)] and 
by definition II contains the necessary countcrtcrm to 
vanish at the true mass k 2 = m 2 . To do perturba- 
tion theory with mass mo the full propagator is written 
D'(k) = l/[fc 2 -m 2 -n(fc)] where fl(k) = m 2 -m 2 +U(k). 
Of course II does not vanish at k 2 = m 2 or at k 2 = m\ 
and this is the source of the problem. A self-energy in- 
sertion on the internal lines of Fig. 1 gives the two-loop 
contribution shown in Fig. 2. 



FIG. 2. Two-loop self-energy due to one insertion. 
The corresponding self-energy is 

n( 2 >(fc) = l9 2 J ^ [A(p)n(p)A(p)]A(p- k). (3) 



This has a two-particle and a three-particle discontinuity. 
The quantity in square brackets has a double pole at 
p 2 = m§ since Tl(p) does not vanish at p 2 = m\. This 
produces a two-particle discontinuity of the form 



Discn( 2 )(fc) = ^- 
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The correction is is infinite at k 2 = 4m 2 , and large near 
there. This is a signal that the correct branch point is 
not at k 2 = Am 2 . Multiple self-energy insertions produce 
successively higher powers of the inverse square root: 
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This is the beginning of a Taylor series. All the correc- 
tions diverge at the false threshold k 2 = 4m 2 ,. In the 
range Am\ < k 2 < Am 2 each correction is finite but the 
Taylor series diverges. Thus perturbation theory fails 
throughout this region of k 2 . To obtain a convergent se- 
ries it is necessary to work in the range k 2 > Am 2 . In 
this region the Taylor series converges and the sum is the 
full two-particle discontinuity with branch point shifted 
to the physical mass m 2 = m\ + dm 2 : 
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The true two-particle threshold is still a square root 
branch point at k 2 — Am 2 . The breakdown of pertur- 
bation theory is entirely due to a propagator A(fc) = 
l/[fc 2 — Wq] with the the wrong mass m . The break- 



down is easily avoided by using 1/[A: 2 
particle propagator. 



m 2 ] for the free 



B. Breakdown of thermal perturbation theory 

In the preceding example mo was the wrong mass. Fi- 
nite temperature calculations are normally done using 
free propagators with the wrong mass, i.e. the zero- 
temperature mass m. A typical component of the con- 
ventional finite-temperature matrix propagator A a b(k) is 



A ++ (k) 



1 



2iriS(k 2 - m 2 ) 
k 2 — m 2 + ie exp(/3|fco|) — 1 



The two-loop topology shown in Fig. 2 produces a matrix 
self-energy 

ni 2 j (k) ^ig 2 /^ [Aa h ( P )n« ( P )A cd ( P )] A ad (k P ). 

As in Eq. (3) the quantity in square brackets has a dou- 
ble pole at p 2 = m 2 and is the source of the problem. 
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The self-energy has branch cuts at two-particle and three- 
particle thresholds as well as cuts due to absorption of 
one or more particles from the heat bath. The disconti- 
nuities across any of these is given by the cutting rules of 
Kobes and Semenoff |6| . A particular linear combination 
of the H a d gives the retarded self-energy. The formula for 
the discontinuity of the two-loop retarded self-energy is 
given explicitly in Jt],|| in terms of the one-loop n^P . It 
is easily evaluated in the limit fc 2 — > Am 2 with the result 



DiscIT^fc)- 



-ig 2 [1 + 2n(fc /2)]Rer4 1) (fc /2) 
16tt fc 2 (l-4m 2 /fc 2 ) 1/2 



The correction is infinite at k 2 — Am 2 indicating the 
breakdown of finite-temperature perturbation theory in 
this region. This breakdown is an artifact of using prop- 
agators with mass m. 



III. QUASIPARTICLE RESUMMATION FOR <j) 3 
FIELD THEORY 

To cure the breakdown of thermal perturbation the- 
ory described above it is necessary to use a free thermal 
propagator that has the same poles as the exact thermal 
propagator. The exact thermal propagator matrix may 
be expressed as a linear combination of the exact retarded 
and advanced thermal propagators. The exact retarded 
thermal propagator D' R (k) has the following properties: 
(a) It is analytic for Im ko > 0; (b) It has a pole in the 
fourth quadrant at fco = £ (fc); (c) It has a pole in the 
third quadrant at fco = —£*(k) because of the relation 
D' R (ko,k) — \D' R {— fc^ fc)] *• It will be useful to have an 
explicit notation for the real and imaginary parts of £: 



£(k) = Elk) - iT(k)/2 



(E>0;T> 0). 



(5) 
iT/2. 



The fourth-quadrant pole is thus at — £* = —E 
The functions E and T also depend on T,g,m. 

It is somewhat misleading to write the retarded prop- 
agator as D' R (k) — l/[fc 2 — m 2 — n^(fc)] since there is no 
pole at fc 2 = m 2 . It is better to express it as 

D' R {k) 



(ko-^iko+^-nRgpik)' 



where U Rqp (k) = U R (k) + fc 2 + m 2 - \£\ 2 + iTk a . It is 
natural to treat HR qp perturbatively since in vanishes at 
fco = £ and at fco = — £* ■ Consequently the free retarded 
quasiparticle propagator is 



D R (k) 



(k - £)(k Q + £*) 
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The exact advanced propagator is D' A (k) = D'A—k) and 
has poles in the first quadrant at fco = £* and in the 
second quadrant at fco = —£■ Consequently the free ad- 
vanced quasiparticle propagator is 



D A (k) 
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It is straightforward to reorganize the Feynman rules 
for perturbation theory so as to use the free quasiparticle 
propagators. The 2x2 propagator matrix is 



D ab (k) 



D R (k) 
-D A (k) 



(8) 



+[D R (k)-D A (k)}n(k ){ ( ;_ CT)fc0 1 y 

where n(fco) = l/[exp(/3fco) — 1] is the Bose-Einstein func- 
tion without absolute value bars. This satisfies the KMS 
condition jj] . Whenever there is a self-energy insertion on 
a propagator (i.e. DUD) then a counterterm <5II should 
be inserted where 



SU ab (k) = (k 2 +m 2 -£*£) 



1 

-1 



(9) 
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The check of this counterterm is as follows. The inverse 
of Eq. (8) is 

D- b \k) = {kl-£*£)(l -1 

/ e P k ° -f i — 2e ak ° 

+ik Q Tn(k ) f _2 e (/3-er)fc e f3k„ _|_ ]) " 

(10) 

The identity 



D- b \k)-m ab (k)=( *V 2 _/ +m2 



guarantees that the two approaches are identical when 
summed to all orders. 



A. One-loop retarded self-energy 

Using the quasiparticle propagator (||) it is straightfor- 
ward to compute the one-loop self-energy in Fig. 1. The 
retarded self-energy reduces to 



n*(fc) = \ 



g 2 fd 3 k 1 d 3 k 2 



S 3 (h+k 2 -k)/(2TT) 3 (11) 



2 J 2E 1 2E 2 

/[l + n(£i)][l + n(g 2 )]-n(giMg2) 



fco - £\ - £2 




[l + n(£i)]n(££)-n(£i)[l + 


ni£l)\ 


fc - £1 + £ 2 




n(£f)[l + n(£ 2 )]-[l + n(£f[ 


)]n{£ 2 ) 


fco + £{ - £ 2 




n(£t)n(£i)-[l + n(£ ?)][! + 
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fco + £{ + £* 2 
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The statistical factors in the numerators of Eq. 
describe the emission and absorption of quasiparticles, 
weighted by 1 + n(£) and n{£) respectively, and also 
the emission and absorption of quasiholes, weighted by 
1 + n(£*) and n(£ *). Each of the numerators describes 
a direct process minus the inverse process. The denom- 
inators are singular in the lower half-plane at the four 
locations 

k = X(h) + X{k 2 ) X = £or-£*. (12) 

All of these lie in the lower half-plane since £ and — £* 
have the same negative imaginary part. (There are also 
unphysical cuts related to Matsubara frequencies but 
these are cancelled by the two-loop contribution Eq. 
(11) can also be written as a dispersion integral over a 
spectral function as discussed in Eq. (3.124) of [10]. 

B. Two- loop retarded self-energy 

The two-loop self-energy shown in Fig. 2 requires the 
counterterm matrix. The retarded part part of the self- 
energy has both two-particle and three-particle cuts. The 
two-particle cuts are not singular (because the integrands 
contain no double poles) and have the same form as Eq. 
( |TT| ) only multiplied by wave-function renormalization 
factors. There is no breakdown of perturbation theory 
near the production thresholds because the location of 
the thresholds is now correct. The three particle cuts of 
Fig. 2 occur at 

k Q = X(k 2 ) + A(fc 3 ) + X(k 4 ). (13) 

There are eight possibilities since A — £ or — £*, inde- 
pendently. 

The two-loop retarded self-energy shown in Fig. 3 has 
a variety of branch cuts. It has two-particle cuts at 

k = Mh) + x(k 2 ) 

k a = A(fc 3 ) + A(fc 4 ), (14) 

and three-particle cuts at 

k = X(h) + X(k 3 ) + X(k 5 ) 

k =X(k 2 ) + X(k 4 ) + X(k 5 ). (15) 

A complete account of the two-loop calculations is given 
in|. 
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Fig 3: Two-loop self-energy due to vertex correction. 



IV. FIELD THEORY OF QUASIPARTICLES 

Reorganizing the Feynman rules so as to use the free 
quasiparticle propagator (||) accomplishes two things. 
First, the perturbative expansion does not break down 
near the branch points. Second, the branch points of in- 
dividual diagrams are the correct branch points of the 
full theory. The original motivations given in Sec. I can 
now be restate more specifically: 

• Do the quasiparticle Feynman rules come from an 
operator theory with a Hamiltonian that acts in a 
linear vector space? 

• Are the poles in the retarded propagator at ko = 
£(k) and at ko = —£(k)* eigenvalues of a Hamilto- 
nian? 

• Are the energies of two-particle branch cuts at 
fc = X{kx) + X{k 2 ), where A(fc) = £{k) or -£*{K), 
eigenvalues of states containing two quasiparticles? 

The answer to all these questions is "yes" . However since 
the energies A(fc) are temperature-dependent and com- 
plex the relevant Hamiltonian cannot be the usual one. 

Previous treatments have not used the complex energy 
pole to define quasiparticles, but instead have used a con- 
tinuous mass spectrum P~pT| . These approaches are not 
closely related to methods for calculation as it is unclear 
how to incorporate interactions. 

The treatment of quasiparticles described below is 
based on a Lagrangian L that enjoys two properties. 
First, it generates the conventional finite-temperature 
perturbation series with free propagator poles at k 2 = 
to 2 . Second, by adding and subtracting the same terms 
the Lagrangian is written as L = Lq + Li in which Lq de- 
scribes free quasiparticles with the propagator (|J) and Lj 
generates the counterterm SH and the self interactions. 

A. Lagrangian for free quasiparticles 

The conjecture for the free quasiparticle Lagrangian 
is based on the fact that the matrix propagator (g) has 
poles at four distinct locations: £,—£*,—£,£*. A field 
equation that is second order in the time derivative can 
only contain two of these energies. If the fields are self- 
adjoint, the only possible self-adjoint equations are 

(i- + £*)(-i-+£)ct> A (t,k)=0. (16) 
The Lagrangian that produces these equations is 

^ [H^+n<M«, -*)] [(< JpOM*, *o] ■ 
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To convert the integration over k to an integration 
over spatial coordinates requires introducing the spatial 
Fourier transform of the width T(k) given in (5), 



T{x-x" ) 



d 3 k 
(2tt) 3 



e ifc. r(fc) 



(17) 



and similarly for £*£(k). The Lagrangian becomes 

L a(t)= j d 3 x 4> A (x)(f> R (x) 

+ J d 3 xd 3 x' [- <f> A {x)£*£(x - x')<f> R (x') (18) 

~$ A {x)T(x - x')(f> R (x') + ~(/> A (x)T(x - f )j> R (x')]. 

The fields are all evaluated at time t = xq = x' . 

Before examining how to quantize Lqq, it is helpful to 
present an alternative argument that provides the rela- 
tion between the fields <f> R / A and the original fields. The 
free quasiparticlc action may be written as 



dt Lqo — 



d 4 k 



M-k)(ko+£*)(ko-8)<l>A(k). 



This form is closely related to the results of Aurenche 
and Bccherrawy jjjj and of van Eijck, Kobes, and van 
Weert ]l3|,0| on transforming the matrix propagator (||) 
to a skew-diagonal form. In the notation of JLJ] the prop- 
agator may be 'diagonalized' as 

D ab (k) = VT(-k )( D ° {k) Da W^ V 7b (k ). (19) 

Latin letters a, 6, c . . . denote the original ± components; 
Greek letters a,/?, 7,... denote R/A. The matrix 
has diagonal entries D RR — D AA = and off-diagonal 
entries D RA = D A , D AR = D R . The matrix Vph(ko) has 
rows labelled by R/A and columns labelled by ±. The 
inverse of the center matrix D^(k) in Eq. (19) is 



(k -£)(k + £*) 

(k -£*){k + £) 



This is a Hermitian matrix. Using the doublet notation 



the free quasiparticle action can be expressed as 



dt Lqq — — 



d A k 



jM-^D-^Mk). (20) 



This formulation of the action suggests that the usual 
fields 4>b (b = ±) be related to the new fields 4> 1 (7 = 
R/A) by the matrix transformation 



4> b (k) ~ 4>-y b (k)Vy(k ). 



The transformation is non-local in time, 

/oo 
dt'(j> 7 (t +t',x)V yb (t'). (22) 
-00 

and is certainly not a canonical transformation. The 
symbol ~ does not imply mathematical equality. How- 
ever it will suggest the full Lagrangian (26) below. 



B. Full Lagrangian 

The full Lagrangian will determine the evolution of 
the fields 4> R j A . It is helpful to examine the contour La- 
grangian used in the conventional path integral formula- 
tion of real-time, finite-temperature field theory: 



1,2/2 9 ,3 

2 ^ ~ 3! 0+ 



-\{d^-)(d^_) + ^m 2 4>l + 

Under the transformation (22) the quadratic terms of C c 
become 

{d^ R ){d^ A )~m 2 M A . 



However under (22) the cubic powers (4>±) 3 contain prod- 
ucts of the fields <fi R / A at different times. To avoid this it 
is necessary to make a truncated Taylor series expansion 



JV-l 



d^yix) . 
dxf, 7± ' 



where the coefficients are temporal moments: 



a 



>, - \ dt' -V lb {t')/T{l + ai). 



(23) 



(24) 



(The parameter a guarantees convergence later. It has 
positive real part. At the end of any calculation a — > 0.) 
With the definition (23), cubic powers of <f>± will only in- 
volve products of fields evaluated at the same time. Cal- 
culations must be done at arbitrary N and subsequently 
the limit N— > 00 is performed. In this limit the moments 
have the the property 



lim > (ikn) 
a->o ^ 



V lb {k ) 



(25) 



1=0 



which will be useful later. 

When the definition (23) is substituted into C c the 
Lagrangian will contain time derivatives up to order N-l. 
Consequently there will be N-l canonical momenta. To 
construct the Hamiltonian requires inverting the relation 
between the highest time derivative of the field and the 
highest momentum. This inversion is not tractable if the 
highest time derivative is buried in the series definition of 
(21) <j)±. The solution is to add a quadratic term containing 
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a time derivative of order N. The full Lagrangian density 
is then 



m 2 (fi R (f> A 

M 3 . 



d 



N, 



! >R 



d 



N, 



dxtf 



N 



(26) 



Because of the auxiliary term it is trivial to invert the re- 
lation between the highest time derivative and the N'th 
canonical momentum. In order that h has no physical ef- 
fects it will have to approach zero as N^oo. The specific 
choice is 



h N 



-JV 



(27) 



V. QUANTIZATION OF Zoo 

Although the unperturbed Lagrangian L$ produces 
field equations with time derivatives of order 2N, the 
vanishing of h according to Eq. (27) is chosen so that 
as N— > oo the extra modes become infinitely heavy and 
decouple ||. Consequently the essential physics of the 
free quasiparticles is contained in Loo- It it therefore 
quite instructive to quantize Loo first. This will lead to 
a self adjoint Hamiltonian Hqq whose eigenvalues are the 
locations of the poles and cuts in the exact thermal prop- 
agators. 



A. Field operators 



C. Comments on quantization of L 

The Lagrangian density (26) has all time derivatives 
up through order N. The construction of a Hamiltonian 
for such systems was developed by the mathematician 
Ostrogradski in 1850 and is discussed in the mechanics 
book by E.T. Whittaker @. For the full Lagrangian 
density the 2N canonical coordinates are the fields 4> r /a 
and their time derivatives up through order N-l. There 
are 2N canonical momenta. The classical Poisson bracket 
is replaced by the quantum mechanical commutator in 
order to quantize. 

To compute in perturbation theory it is necessary to 
choose the unperturbed Lagrangian. It must include the 
term proportional to h in order for all the canonical mo- 
menta to appear in the unperturbed theory. One choice 
for the unperturbed Lagrangian density would be the 
three quadratic terms on the first line of Eq. (26). This 
choice ultimately leads, in the limit N— > oo, to conven- 
tional Feynman rules with propagators containing poles 
at k 2 = m 2 as shown in || . 

For the unperturbed operator theory to describe quasi- 
particles, a different choice of unperturbed Lagrangian is 
necessary: 



Lq — Loo — h 



d N R d^A 

dx$ dxg 



(28) 



This also has 2N canonical coordinates and 2N canonical 
momenta. In the limit N— > oo the free propagators will 
be D R and Da- The full Lagrangian (26) is separated 
into L = Lq + Lj, where the interaction term is 

L,=J d 3 z[- ■ V0 fl - m 2 A 0R - j { 4>1 + 



d 3 xd 3 x' 



(t> A {x)e*£{x-x H ) ( t> R {x') 



(29) 



+ -4>a{x)Y{x - x)4>r{x) - -4>a(x)T(x - x)4> R {x) 

The quadratic terms in Lj will produce the counter term 
matrix (9) and the cubic terms produce the interaction 
vertices. 



The field equations that follow from Loo are given in 
Eq. (16). The solutions can be expanded in terms of 
spatial Fourier transforms in a box of volume V with 
periodic boundary conditions: 



<p R (t,g) =22 



^ t R ik-x — i£*t | R\ —ik-x J ri£t\ 
k 



(30) 



(f> A (t,x) 



1 



(qA ^ik-x — i£t j_ ^A~\ ^ — ik-x-\-i£*t\ 

V2EV * 1 ' 

k 



Note that <f> R (t, fc)<~ exp(+rt/2) rises exponentially with 
time and ^(t, exp(— Tt/2) falls exponentially. The 
Hamiltonian will contain products of the two fields and 
will be time-independent. For a fixed wave vector fc, one 
may interpret <pA{t, k) and 4> R (t, k) as the coordinates of 
a damped harmonic oscillator an anti-damped harmonic 
oscillator, respectively. The oscillator problem, devoid of 
the field concept, has a long history. Bateman |L6| in 
1931 invented a Lagrangian of the form (18) to treat the 
classical problem of one damped and one anti-damped 
oscillator. Morse and Feshbach |l7j also use this same 
oscillator Lagrangian. The quantum mechanics of the 
system is discussed in M,19l . 



B. Canonical momenta and commutation relations 

Since there are only first order time derivatives, the 
usual canonical formalism applies directly. The conjugate 
momenta are 7r Q = SLoo/54> a '- 

ir R (x) = 4> A {x) + d 3 x'T(x-x')<j)A{x') 

n A (x) = <p R {x) - i J d 3 x'T(x - x')(p R {x'). 

The canonical commutation relations are 



= -S a0 S 3 (x-y). (31) 

xo=yo i 
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Consequently the raising and lowering operators satisfy 



At 



1 fc 





and 



The other commutators vanish: 
and [a R ,a A ] = 0. This implies [4> R (x),4> R (y)} 
[<f>A{x),<t>A{y)] = 0. Because the time dependence of (f> R 
and (j>A, if these two field commutators did not vanish 
they would have to depend on both xq and yo (rather 
than x — ?/o ) an d would thus violate invariance under 
time translations. The only non- vanishing commutator of 
two fields is [<j>A{%)i<j>R,{y)] an d it automatically satisfies 
translational invariance. 



C. Hamiltonian Hoo 



Note that both £ and — £* have the same negative imagi- 
nary part. The next excitations are those with two quasi- 



(32) particles 



a 1 a 



\0r), with one quasiparticle and one 



quasiholc, a R ^a R \0 R ), or with two quasiholes, a R a R \0 R ). 
Every state of the form (a R ^ ) e (a R )' m |0 R ) is of the re- 
tarded type and has an energy with a negative imaginary 
part. (The real part of the energy may be of either sign.) 

There is another set of states that are excitations of 
the advanced vacuum state \Qa), which is defined by the 
conditions 



of|0x> 







af 10.4) =0. 



(35) 



The advanced representation of the one quasiparticle ex- 
citation is a^lOyi) and of the one quasiholc excitation, 
cl a \0a)- From the commutation relations (32) these op- 
erators shift the energy by £* or —£, respectively: 



-k r $ r + tt a 4>a] 



00 



The Hamiltonian iJoo = J d 3 x 

describes free quasiparticles. In terms of canonical vari- 
ables it is 

H a = J d 3 xd 3 x'[n A (x)S(x - x)-k r {x) 
-]^a{x)Y{x - x')(f> A (x') 
+ ^ r {x)T(x-x")k r {xI) 
+Mx)E 2 (x-x')4> a (x') 



The kernel E 2 (x — x') is the Fourier transform of E 2 (k) 
analogous to (17). When the plane wave expansions (30) 
are substituted the Hamiltonian reduces to 



, {) =^[£ {a fa A + -)+£* 



Hnn — 



(33) 



which is manifestly self-adjoint. The complete analysis 
of this Hamiltonian is rather lengthly. The results of 
that analysis will be summarized here but not proven. 
The essential feature is that there are two parallel sets of 
eigenstates. One set is built upon a ground state |0.r); 
the other, on a ground state \0a}- The retarded vacuum 
state \0 R ) is defined by the conditions 



af\0 R ) 







af\0 R )=0. 



(34) 



(Recall that a A and a A ^ commute.) The retarded vac- 
uum is an eigenstate of ffoo- (Since H no is not normal 
ordered, the eigenvalue is not zero.) The simplest excita- 
tions of the vacuum are the addition of one quasiparticle, 
a R ^\0 R ), or the addition of one quasihole, a R \0 R ). From 
the commutation relations (32) these operators shift the 
energy by £ or — £*, respectively: 



[ffoo, af] = £af [H 00 



a R ] = -£*al 
fc j fc 



-£a A . 



Both £* and —£ have the same positive imaginary part. 
The advanced representation of a general excitation is 
(a A ^) e (a A ) m \0A) and has an energy with a positive imag- 
inary part. 

As explained in Appendix A, the Hamiltonian docs 
not act in a Hilbert space. None of the retarded states 
\^r) have finite norm; none of the advanced states \^a) 
have finite norm. However the scalar product between 
retarded and advanced states, is always finite. 

For example, (0 A \a A a Rf \0 R ) = 1 and (0 A \a^af\0 R ) = 
— I. Physical applications always involve matrix elements 
of this form. The completeness relation, 



(36) 



is useful in evaluating matrix elements of the form 
(0a|A^F|0,r) as shown in the next section. 



VI. SAMPLE RESULTS 

It is not possible in the allotted space to give a full 
discussion of this formalism. One of the important fea- 
tures that will be used below is that all the temperature- 
dependence is contained in the full Lagrangian (26). Con- 
sequently propagators, Green functions, and the free en- 
ergy are computed as vacuum matrix elements; there are 
no thermal traces. 



A. Free propagators for 4>r/a 

The quasiparticle pole in the propagator comes directly 
from the eigenstate of the Hamiltonian i?oo- The free 
retarded propagator is defined as a time-ordered product: 
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D AR {x) = -i{0 A \T(<f> A (^n(0))\0 R ) 

= -ie(x ){o A \<p A (x)Mo)\OR) 



N-l 



(37) 



Only one time-ordering survives on the second line be- 
cause (0 A \4> R (0)4> A (x)\0 R ) = 0. When the plane wave 
solutions (30) are inserted, the result is 



D AR (x) = -i8(t) Y,^{ e ^ m ~ e-* 5+l£ * 4 



The four-dimensional Fourier transform is 

1 



D AR (k) 



(fc - £)(k + £*) 



(38) 



Thus the operator theory produces the quasiparticle 
propagator as desired. Similarly, the free advanced quasi- 
particle propagator is 



D RA {x) = -i(0 a \T(Mx)M0))\0r) 

= -i6(-x q )(0 a \M0)Mx)\0r) 



(39) 



Comparison with (37) shows that D RA (x) = D AR (—x). 
Thus D RA {k) = D RA (—k). Since (0 A \<b R q> R \0 R ) = 
and (0 A \(j) A (j) A \0 R } = the matrix propagator in the re- 
tarded/advanced basis is 



D a p(k) 



D RA {k) 
D AR (k) 



(40) 



in agreement with Eq. (19). 



B. Free propagators for <f>± 

Although the complex quasiparticle energy £(k) is 
temperature-dependent, there are no Bose-Einstein func- 
tions in the retarded or advanced propagators (40). How- 
ever, the Bose-Einstein functions should appear in the 
propagators for <fi± in order to agree with Eq. (|§|). The 
(j>± propagators are vacuum matrix elements: 



ab (x -y) = ~z(0 A \T^ a (x)My)) |0h), (41) 



G. 



where a, b = ±. To compute this properly requires quan- 
tizing not Loo but Lq given in Eq. (28). I will provide the 
essential feature of that quantization at the appropriate 
point below. 

Using the definition (23) of the <fi± as a sum of the first 
N-l time derivatives of <j) R and 4> A gives for G ab 



G ab (x-y) 



N-l 



In \rrf d M x ) ^ i/ 'im \ n i n m 



When 4> R / A are quantized according to Lq the time 
derivatives may be taken outside of the time-ordering so 
that 



Gab(x - y) = ^2 



d d 



Lmto d < dy r o 



G/3~{(x — y)C f3a C™ b 



The Fourier transform is 

N-l N-l 



m=0 



As the number of derivatives goes to infinity the identity 
(25) allows the moments to be summed: 

^G ab (k) = V a T (-k o ) ( DA ° Rik) D «f^ ) Vfiiho). 



This is exactly the result (|19|) for the matrix propagator 
of quasiparticles. The Bose-Einstein functions are con- 
tained in V(ko). 



C. One-loop self-energy 

There is a pleasant surprise when we examine the inter- 
acting propagator. The one-loop self-energy, for example, 
is 



n bc (a; - y) 



-ig 



a b a c (Q A \T{<t>l{x)4> 2 c (y))\0 R ). (42) 



One way to compute this is to use Wick's theorem to 
express the matrix element as [D bc (x — y)] . Comput- 
ing H bc (k) then requires integrating the one- loop diagram 
Fig. 1 with quasiparticle propagators. The result is al- 
ready given in Eq. (11). 

The operator expression (42) can be evaluated in an- 
other manner that explains all the branch cuts displayed 
in Eq. (11). The quasiparticle Hamiltonian determines 
the time-dependence of the field operators: 



cj> b (x) 



e lHooXo (f> b (x)e- lHooXo 



as can be easily verified. For either of the two time- 
orderings in Eq. (42) a complete set of states may be 
inserted between the operators using the completeness 
relation (36). The states are energy eigenstates satisfying 
HooI^r) = \*r)£$ and H o\^a) = with £$ = 

£§* . The matrix element in Eq. (42) has the value 



£ [6(x - y„) <0a |$ (x) | (0r A \<£ (y)\0 R )e-^»-yo) 
+%o - x ) (0a I 4>l (y) I *b) (*a 1 4 (2) |0 B ) e l£ S («<>-*>) 



Since Im £ R < the time dependence is exponentially 
falling in both terms. The temporal Fourier transform to 
kg at fixed x and y is 



(0 a \<%W\* r )(* a \cI?M\0r) 

<p • fc o - £$ 

(0 A \tf(y)\* R )(* A \rt(x)\0 R ) 



£^ 



(43) 
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The operators , only connect the vacuum with states 
containing two excitations (i.e. two quasiparticles, one 
quasiparticle and one quasihole, or two quasihole). Thus 
the possibilities are 

a^afjOn) £{k x )+£(k 2 ) 

a^afJOn) ( 44 ) 

The four terms in IIf, c with denominators k — £^ all have 
branch cuts in the lower-half of the complex ko plane. 
These produce the four branch cuts previously found in 
Eq. (11). The retarded self-energy is the particular lin- 
ear combination of the Tlb c which contains only these 
four cuts. The four terms in Eq. (43) with denomina- 
tors fco + £ R all have branch cuts in the upper-half of the 
complex ko plane and produce the advanced self-energy. 
The numerators in Eq. (43) may be evaluated using rela- 
tions (23) and (25) and result in the correct Bose-Einstein 
functions n(£) and n(£ *) in agreement with Eq. (11). 



Bose-Einstein function with the T=0 energy as argu- 
ment: n(Eo). Consequently the transformation from the 
original fields (j>± to the TFD fields </> and <j> depends on 
space rather than on time. The transformation is auto- 
matically canonical. See Appendix B of [5]. (When TFD 
is extended to nonequilibrium situations, the transforma- 
tion is time-dependent |l0| , ^3|| .) 

In both formulations the number of degrees of freedom 
is doubled. In thermofield dynamics the operators act 
in a Hilbert space. The basis states are built out of two 
commuting raising operators, and a^. In the approach 
described here, operators do not act in a Hilbert space, 
as discussed in Appendix A. Furthermore, there are two 
parallel sets of basis states \^r) and \^a)- The former 
are built from commuting raising operators a R<! and a R 
acting on |0/j); the latter from a Aj; and a A acting on |0,4). 
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VII. DISCUSSION 

The treatment presented here is based on the fact 
that field theories at finite temperature inevitably have 
poles in the full propagator at complex, temperature- 
dependent energies £{k). To consistently include quasi- 
particles in Feynman diagram computations it is only 
necessary to use the quasiparticle propagator (8) and 
the counter-term (9). To go further and incorporate the 
quasiparticle effects into the operator structure requires 
the temperature-dependent Lagrangian (26). With this 
Lagrangian the questions posed at the beginning of Sec. 
IV may all be answered affirmatively. 

The logical necessity of this development is rather 
tight. As shown in Sec. II, higher loop corrections will re- 
sult in a breakdown of perturbation theory unless the free 
propagators are chosen to have the correct quasiparticle 
poles. The free quasiparticle propagator shown in Eq. 
(8) must then contain the off-shell Bose-Einstein func- 
tion n(fco) in order to satisfy the KMS condition. (With 
quasiparticle poles, neither the quasiparticle energy £{k) 
nor the T=0 energy E Q (k) — (k 2 + m 2 ) 1 / 2 can be the 
argument of the Bose-Einstein function.) The off-shell 
n(ko) leads to the time-dependent relation (22), or the 
more precise relation (23), between the original fields and 
4'r/a- That fixes the Lagrangian to be (26). 

The standard operator formulation of finite tempera- 
ture field theory is thermofield dynamics [p0|-p2[. The 
unperturbed propagators in thermofield dynamics have 



APPENDIX A: COMPLEX EIGENVALUES 
FROM A SELF-ADJOINT HAMILTONIAN 

States in a Hilbert space have finite norms and conse- 
quently self-adjoint operators cannot have complex eigen- 
values. The state space for quasiparticles is not a Hilbert 
space. A retarded cigenstatc of the quasiparticle Hamil- 
tonian satisfies 



HqoI^r) 



\*r)£ r - 



(Al) 



Since Hqq is self-adjoint, the adjoint of this equation is 
(^> R \H ao = £ R *($f R \. The inner product of this with 
|*_r) would force £ R to be real If the norm (^rI^r) 
existed. As shown below, this norm does not exist. 
Advanced eigenstates satisfy Hqo\&a) = \^a)£ A - The 
scalar product of this with r\ does exist and implies 



<S> R \$ A )£ R * = (*r\$a)£ j 



(A2) 



Thus if £ R * 7^ £ A then the states are orthogonal. 

A concrete realization of the quasiparticle algebra is 
illuminating. The four operators a A ,a A \a R ,a R ^ may 
be represented in a two-dimensional real space as linear 
combinations of x, d/dx, y, d/dy as follows 



a A + <r 4t = x - iS- a A - a At = — 



poles at k 



The matrix propagators contain the 



a R + a RJ< = x + i 



dy 

d_ 

dy 



dx 



W 



a R - a fit = f- x - iy (A3) 







These satisfy the commutation relations (32). The free 
quasiparticle Hamiltonian (33) in this representation is 



E , d 2 2 d 2 ~ r, d 2 s 
Hoa = j(-d^ + x + W~ v) + 2 { d^-y +xy) - 



The retarded vacuum state, defined by the conditions 
ftj4 |0fl) = a^|0fl) = 0, has the representation 



(x\0r) = cxp(-2.Ty)/V7r, 



(A4) 



which is not square integrable. The wave function is an 
eigenstate of Hqq with eigenvalue —iT/2. Application of 
the quasiparticle creation operator and the quasihole 
creation operator a R give excited state wave functions 



(x\(a R ^ e (a R ) m \0 R ) = H e (x + iy) H m (x - iy) ■ 



-ixy 



where Hi and H m are Hermitc polynomials. None of 
these are square integrable. These are eigenstates of Hqq 
with eigenvalue £(£+1/2)- £*(m + 1/2). 

The advanced vacuum state is defined by o, r \0a) = 
o> r ^\0a) — and has the representation 



(x\0 A ) = exp(ixy)/^. 



(A5) 



It is not square integrable. However the scalar product 
of the advanced and retarded vacuum states is finite: 



(0a\0r) 



dxdy ex PH 2 ^) = L 



(A6) 



The excited state wave functions are 



(x\(a^Y(a A ) m \0 A ) = H e (x - iy) H m (x + iy) 
All scalar products (^a\^r) are finite. 
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